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ABSTRACT 


The  purpose  of  this  research  was  to  compare  numerical  predictions  with  experimental 
results  and  to  devise  an  accurate  laser  anemometiy  technique  to  measure  in  the  endwaU 
region  of  a  confined  annulus.  Flowfield  characteristics  were  computed  using  a  three- 
dimensional  flow  solver  with  the  numerical  plane  coincident  with  the  experimental 
measurement  location.  An  annular  turbine  cascade,  designed  for  laser-Doppler  velocimetry, 
was  modified  to  obtain  blade  passage  midspan  surface  pressure  measurements.  A  range  of 
predicted  subsonic  and  transonic  midspan  surface  pressures  were  compared  favorably  with 
experimental  measurements.  Two-dimensional  Mach  number,  flow  angle,  and  turbulence 
intensity  measurements  were  obtained  with  a  fiber-optics  laser-Doppler  velocimeter.  The 
measurements  were  performed  through  a  1 .0922  millimeter  opening  in  the  endwall  at  depths 
ranging  from  0.01  mm  to  3.34  mm  and  the  results  were  compared  with  numerical 
predictions. 
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I.  INTRODUCTION 


Advanced  aeropropulsion  systems  have  required  substantial  technological 
improvements  in  turbomachinery.  These  depend  on  the  availability  of  complex  analysis  and 
experimental  tools.  A  thorough  understanding  of  the  generation  and  development  of 
secondary  flows  in  annular  blade  rows  has  provided  insight  into  the  design  and  performance 
of  turbomachinery.  Secondary  flows,  together  with  tip  leakage  flows,  produce  considerable 
flow  distortions  and  losses  in  the  endwall  region.  [Ref.  1]  These  losses  can  be  minimized  and 
turbomachinery  efficiency  can  be  improved  with  a  more  precise  understanding  of  the  flow 
mechanics  and  the  ability  to  numerically  predict  the  flow  field.  The  current  emphasis  on 
turbomachinery  design  centers  on  numerical  analysis. 

This  report  specifically  documents  the  investigation  of  the  flow  through  an  annular 
turbine  cascade  (ATC).  References  2  through  4  include  continuing  research  in  the  field  of 
laser-Doppler  velocimetry  (LDV)  and  numerical  prediction  in  a  large  annular  turbine 
cascade.  The  authors  provided  experimental  two  and  three-dimensional  velocity,  flow  angle, 
and  turbulence  intensity  at  constant  axial  positions  within  an  ATC  passage.  These 
measurements,  in  addition  to  vane  surface  static  pressure  measurements,  supplied  a  test  case 
for  three-dimensional  turbomachinery  computer  programs.  They  concentrated  on  obtaining 
comparisons  ahead  of,  inside,  and  downstream  of  the  blade  passage  and  concluded  that  the 
largest  difference  between  experimental  and  computational  results  was  in  the  endwall  region 
where  viscous  and  secondary  flow  effects  were  the  greatest. 

Two  previous  investigations  of  the  flow  through  the  ATC  have  been  conducted  at  the 
Naval  Postgraduate  School.  Reference  5  included  design  and  manufacturing  information  of 
the  annular  turbine  cascade  that  was  developed  to  determine  the  limitations  of  LDV 
measurements  in  a  confined  annulus.  Reference  6  included  additional  laser  and  pressure 
probe  access  modifications  and  initial  LDV  measurements  to  the  same  ATC.  This  report 
includes  further  ATC  modifications  for  midspan  blade  surface  pressure  measurements  and 
LDV  measurement  techniques  of  the  endwall  flow.  Radial  two-dimensional  fiber-optic  probe 
traverses  were  performed,  through  a  small  access  hole  in  the  outer  casing,  to  coincide  with  the 
numerical  exit  plane.  Circumferential  surveys  were  obtained  at  different  radial  locations  close 
to  the  endwall.  Blade  midspan  surface  pressures  were  measured  within  one  blade  passage  at 
various  inlet  total-to-downstream  hub-static  pressure  ratios.  Blade  surface  pressure  and 
endwall  flow  measurements  were  compared  with  numerical  predictions  obtained  using  a 
three-dimensional  viscous  computer  program. 
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The  two  most  noteworthy  comparisons  were  at  a  subsonic  pressure  ratio  and  a 
pressure  ratio  corresponding  to  sonic  exit  conditions.  Comparisons  with  the  LDV  data  were 
perforaied  at  the  subsonic  flow  condition.  The  numerical  blade  surface  pressure  distributions 
compared  well  with  the  experimental  results,  particularly  for  the  sonic  exit  condition  for 
which  trailing  edge  shocks  were  predicted. 
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n.  EXPERIMENTAL  SETUP 


Figure  1.  Front  View  of  Experimental  Apparatus 

Flow  stagnation  pressure  was  measured  at  two  upstream  locations.  One  combination 
probe  provided  pressure-setting  information  to  a  mercury  manometer  board  and  a  digital 
readout  of  flow  stagnation  temperature,  while  the  second  probe  was  connected  to  a  Scanivalve 
(Figure  2,  lower  right).  Four  (averaged)  upstream  static  ports  and  four  (averaged)  inner  hub 


A.  TEST  FACILITY  AND  ANNULAR  TURBINE  CASCADE 

Airflow  for  the  annular  turbine  cascade  experiment  was  provided  by  a  VA-312  Allis- 
Chalmers  12-stage  axial-flow  compressor  located  at  the  Tuibopropulsion  Laboratory  of  the 
Naval  Postgraduate  School.  The  compressor  was  operated  at  12,000  rpm  at  various  discharge 
pressures  and  provided  metered  air  to  a  plenum  chamber.  Air  from  the  plenum  was  routed  to 
a  232.918  mm  (9.170  in)  diameter  beUmouth  and  test  section  through  honeycomb  flow 
straighteners  in  a  254  mm  (10  in)  flanged  steel  pipe  as  shown  in  Figures  1  and  2. 
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Figure  2.  Top  View  of  Experimental  Apparatus 

Atomized  glycerin  particles  constituted  the  LDV  seed  material  which  were  introduced 
through  a  7.938  mm  (0.313  in)  diameter  copper  tube  approximately  108  tube  diameters 
upstream  of  the  test  section.  Seed  atomization  was  performed  using  a  commercial  TSI,  Inc., 
Six-Jet  Atomizer  connected  to  the  laboratory  compressed  air  supply.  Laser  traverse 
mechanism  buffet,  caused  by  the  exiting  flow,  was  minimized  with  a  prefabricated  aluminum 
cone  (Figure  2,  upper  center)  and  Reference  6  showed  that  flow  characteristics  were  unaltered 
with  the  cone  attached. 

Each  blade  was  designed  with  a  combination  of  simple  circular  arcs  and  line  segments 
and  included  a  leading  edge  radius  of  2.413  mm  (0.095  in),  trailing  edge  radius  of  0.305 
mm  (0.012  in),  and  axial  chord  of  24.77  mm  (0.975  in).  The  annular  stator  row  was 
manufactured  from  2218-T61  aluminum  and  consisted  of  31  blades  with  a  midspan  spacing 
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downstream  (one-half  axial  chord)  static  ports  were  also  cormected  to  a  Scanivalve  for 
automated  pressure  data  acquisition. 


Figure  3.  Blade  Leading  Edge  View  With  Pressure  Measurement  Modifications 


of  21.77  mm  (0.857  in),  resulting  in  a  blade  solidity  of  .1.14.  The  inner  hub  radius  was 
98.93  mm  (3.895  in)  and  the  outer  case  radius  was  116.46  mm  (4.585  in)  with  the  same 
profile  at  all  radii.  Reference  5  included  the  original  set  of  manufacturing  drawings  and 
Reference  6  included  a  description  of  the  wake  positioning  system. 


B.  PRESSURE  MEASUREMENT  MODIFICATIONS 

The  original  ATC  did  not  provide  the  capability  to  measure  blade  surface  pressures. 
Within  a  single  passage,  seven  suction-side  static  ports  and  four  pressure-side  static  ports  were 
drilled  orthogonal  to  the  blade  surface  at  midspan.  Figures  3  and  4  show  front  and  rear  views 
of  the  instrumented  blades.  As  shown  in  Figure  5,  each  static  port  was  0.406  mm  (0.016  in) 
in  diameter  and  each  spanwise  hole  was  1.321  mm  (0.052  in)  in  diameter. 

All  spanwise  holes  were  sealed  at  the  tips  and  stainless  steel  tubes  were  cemented  into 
the  hub  openings.  The  tubes  were  connected  to  the  Scanivalve  with  plastic  tubing  which  was 
fed  out  through  a  sting  from  the  center  body.  (Figures  1  and  2)  Each  port's  circumferential 
position  was  measured  using  a  2.375  mm  (0.0935  in)  diameter  stylus  and  later  converted 
graphically  to  an  axial  chord  position  to  allow  computational  comparisons.  Figure  5  shows 
port  numbering  and  dimensions. 
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Figure  4.  Blade  Trailing  Edge  View  With  Pressure  Measurement  Modifications 


Figure  5.  Blade  Static  Port  Numbering  Sequence  And  Dimensions 
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C.  PRESSURE  DATA  ACQUISITION 

The  data  acquisition  system,  for  the  pressure  measurements,  is  shown  schematically  in 
Figure  6.  All  data  acquisition  was  remotely  controlled  by  a  Hewlett-Packard  9000  computer 
system.  Appendix  A  contains  the  program  utilized  to  conduct  all  pressure  data  acquisition. 

A  Scanivalve  was  connected  to  a  Model  HG-78K  Scanivalve  controller,  which  in  turn  was 
connected  to  a  Hewlett-Packard,  Model  3456A  Digital  Voltmeter  and  Model  3495A  Scanner 
via  a  HP-IB  instrument  bus.  Scanivalve  calibration  was  performed  to  within  an  accuracy  of 
+/-  0.1  inches  mercury.  Table  A1  in  Appendix  A  relates  each  Scanivalve  port  to  its  respective 
pressure  measurement. 


Figure  6.  Pressure  Data  Acquisition  Schematic 
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D.  LASER-DOPPLER  VELOCIMETER 

Part  of  the  laser  apparatus  is  shown  in  Figure  7.  The  probe,  processor,  computer,  and 
traverse  mechanism  are  shown  in  Figures  1  and  2.  The  LDV  system  is  shown  schematically  in 
Figure  8.  A  LEXEL  Model  95  four- Watt  argon-ion  laser  was  connected  to  a  TSI,  Inc.,  Model 
9201  ColorBurst  multicolor  beam  separator.  The  beam  separator  divided  the  incoming  light 
into  shifted  and  unshifted  beams,  with  the  shifted  beam  receiving  a  40  MHz  frequency  shift 
from  a  Bragg  ceU.  The  two  beams  were  further  split  into  three  polarized  pairs:  green  (514.5 
nm),  blue  (488  nm),  and  violet  (476.5  nm). 

Individual  couplers  on  the  ColorBurst  directed  each  beam  to  the  laser  probe  via  a 
fiber-optic  cable.  Each  fiber-optic  probe  contained  receiving  optics  which  directed  the  return 
signal  to  a  TSI,  Inc.,  Model  9230  ColorLink  multicolor  receiver.  The  ColorLink  provided 
photomultiplier  and  frequency-shifting  functions.  All  conditioned  ColorLink  signals  were 
sent  to  a  TSI,  Inc.,  IFA-750  digital  burst  correlator  where  valid  Doppler  signals  were 
identified  and  digitized. 

The  fiber-optic  probes  were  mounted  to  a  LINTECH,  Model  41583  traverse  table. 

An  Applied  Motion  Products  System  1618  traverse  controller  was  used  manually  to  control 
traverse  table  movement.  All  ColorLink  and  IFA-750  functions  and  LDV  data  processing 
were  accomplished  remotely  by  computer  using  TSI's  menu-driven  software,  FIND  (FLOW 
INFORMATION  DISPLAY)  veision  4.04. 

E.  EXPERIMENTAL  PROCEDURE 

1.  Midspan  Surface  Pressure  Measurements 

Midspan  surface  pressure  measurements  were  obtained  with  the  pressure  data 
acquisition  system.  The  pressure  ratio  (Prat)  was  defined  as  the  downstream  hub-static 
pressure  (Phub)  divided  by  the  upstream  stagnation  pressure  (Pq).  Each  pressure  ratio  was  set 
by  metering  the  upstream  stagnation  pressure  until  a  desired  mercury  manometer  column 
height  was  achieved.  Five  pressure  ratios  (0.5070,  0.6041,  0.6815,  0.8077,  and  0.9054)  were 
considered  and  during  each  run  aU  the  blade  surface  pressures,  Pq,  and  Phub  were  recorded. 

2.  Laser  Alignment 

LDV  alignment  for  endwaU  flow  measurements  was  accomplished  as  shown  in  Figure 
9.  The  objective  of  the  LDV  alignment  procedure  was  to  center  the  probe  volume  in  the 
1.0922  mm  optical  access  hole  at  a  known  and  repeatable  radial  distance.  The  two- 
dimensional  fiber-optic  probe  was  attached  to  a  mounting-bracket  micrometer  which  allowed 
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Figure  7.  LDV  Bread  Board  (with  Laser,  Color  Separator,  and  Receiving  Optics  Module) 


Figure  8.  LDV  System  Schematic 
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probe  radial  travel  in  0.01  mm  increments.  After  zeroing  the  micrometer,  the  traverse  table 
was  manually  advanced  forward  until  the  four  beam  separation  was  minimized,  yet  discernible 
with  the  naked  eye  (Figure  9,  sketch  A).  The  digital  position-indicator  reading  on  the 
traverse  controller  was  noted  and  the  process  repeated  for  forward  travel  imtil  the  same  image 
reappeared.  The  midpoint  of  the  traverse  table  positions  was  defined  as  the  center  of  the 
probe  volume  with  the  face  of  the  optical  access  plate  as  a  radial  reference  point.  Reference  6 
described  complete  dimensions  of  the  optical  access  plate. 


1 

1 

Optical  Access  Hole 

1 

1.0922  mm  diameter 

4  / 

i 

1 

1 

1.14  mm  plate  thickness 

k 

1 
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3 

1 

1 

1 

A 

B 

Figure  9.  LDV  Alignment  Schematic 

Horizontal  and  vertical  alignment  positions  (Figure  9,  sketch  B)  were  obtained  by 
noting  the  digital  position  indicator  on  the  traverse  controller  as  the  probe  volume  touched 
the  left(l),  iight(2),  bottom(3),  and  top(4)  inner  edge  of  the  optical  access  hole.  The 
horizontal  and  vertical  reference  positions  were  defined  as  the  center  of  the  optical  access 
hole.  The  probe  was  then  traversed  inward  by  1.14  mm  (the  thickness  of  the  optical  access 
plate),  at  which  point  the  center  of  the  probe  volume  coincided  with  the  outer  (case)  wall  of 
the  turbine  cascade.  This  alignment  technique  was  repeatable  and  ensured  that  the  probe 
volume  passed  cleanly  through  the  center  of  the  optical  access  hole. 

3.  Endwall  Measurements 

AU  endwall  measurements  were  conducted  at  a  pressure  ratio  of  0.9054.  The  laser 
beams  were  aligned  with  the  downstream  optical  access  hole  as  described  above,  and  then 
traversed  manually  using  the  mounting-bracket  micrometer.  Radial  endwall  surveys  at  0.01, 
0.06,  0.18,  0.42,  0.89,  1.78,  and  3.34  millimeters  from  the  tip  casing  were  conducted  for 
peripheral  (wake)  angular  settings  ranging  from  -8  degrees  to  +8  degrees.  One  degree 
(wake)  position  increments  were  achieved  by  circumferentially  rotating  the  blade  row  and 
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center  body  within  the  outer  casing.  At  each  circumferential  position  the  inner  section  was 
secured  in  place  with  a  locking  bolt  arrangement,  the  design  of  which  was  documented  in 
Reference  6.  Seeding  position,  ColorLink,  and  IFA-750  settings  were  adjusted  for  an 
optimum  LDV  data  rate  and  minimum  noise.  All  settings  ensured  a  minimum  data  rate  of 
100  samples  per  second,  however;  LDV  data  rate  was  extremely  sensitive  to  position  of  the 
wand  which  introduced  seed  into  the  flow. 

A  'random*  mode  processor  setting  allowed  a  total  of  1,024  samples  between  the  green 
and  blue  channels  with  no  user  control  over  sample  distribution.  The  'coincidence'  mode 
acquired  1,024  samples  for  each  channel.  Repeatability  measurements  for  endwall  velocity, 
flow  angle,  and  turbulence  intensity  were  separately  conducted  at  three  wake  positions  for 
both  modes.  A  minimum  pressure  ratio  was  determined  before  data  rate  conditions  became 
unacceptable.  Raw  data  were  converted  within  FIND  and  manually  transferred  to  a 
spreadsheet  for  further  processing. 
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m.  COMPUTATIONAL  FLUID  DYNAMICS 


A.  GRID  GENERATION 

Grid  generation  was  completed  using  the  FORTRAN  language  program 
"Turbomachinery  C  GRID  (TCGRID)"  [Ref.  7].  TCGRID  inputs  consisted  of  four  lines  of 
namelist  inputs  followed  by  a  title,  hub  and  tip  geometry,  and  blade  geometry.  The  blade 
geometiy  was  input  in  cylindrical  coordinates  (z,  theta,  r)  starting  at  the  blade  trailing  edge 
and  wrapping  clockwise.  The  blade  inputs  were  completed  in  stacked  sections  from  the  hub 
to  the  tip.  The  grid  used,  with  resolution  of  150  x  31  x  65  (i,  j,  k),  was  the  same  as  that 
generated  in  Reference  6.  The  i-index  was  defined  clockwise  from  the  lower  (pressure 
surface)  exit  to  the  upper  (suction  surface)  exit,  the  j-index  was  defined  from  the  blade 
surface  to  the  periodic  boundary,  and  the  k-index  was  defined  from  the  hub  to  the  tip.  All 
computational  solutions  were  based  on  this  gird.  The  final  grid  is  shown  in  Figure  10. 
Appendix  B  contains  the  grid  namelist  input  file  and  Figure  B1  shows  the  blade  geometry. 

B.  COMPUTATIONAL  SCHEME 

Flowfield  solutions  were  obtained  using  two  versions  (206  and  208)  of  "Rotor  Viscous 
Code  3-D  (RVC3D)",  a  FORTRAN  language  program  designed  for  analysis  of  three- 
dimensional  viscous  flows  in  turbomachinery.  RVC3D  was  written  to  solve  the  thin-layer 
Navier-Stokes  equations  in  Cartesian  coordinates.  The  equations  were  discretized  using 
second-order  finite-differences  in  space  and  solved  in  time  with  a  fourth-order  Runge-Kutta 
scheme.  Streamwise  viscous  terms  were  neglected  using  a  thin-layer  assiunption,  but  cross¬ 
channel  viscous  terms  were  retained.  A  spatially-varying  time  step  and  implicit  residual 
smoothing  were  used  to  accelerate  convergence.  [Refs.  8-10]  Turbulence  effects  were 
modeled  using  a  3-D  adaptation  of  the  Baldwin-Lomax  model  and  the  Cebeci-Smith  model. 

RVC3D  version  206  was  used  to  predict  the  flowfield  for  pressure  ratios  of  0.6041, 
0.6815,  0.8077,  and  0.9054.  The  turbulence  in  the  flowfield  was  computed  with  an 
adaptation  of  the  Cebeci-Smith  turbulence  model.  This  version  of  RVC3D  only 
accommodated  subsonic  exit  boundary  conditions  whereby  the  hub  static  pressure  was  held  at 
the  pressure  ratio  and  radial  equilibrium  was  solved  for  the  spanwise  pressure  distribution. 
Version  208  of  RVC3D  was  used  for  the  0.5070  pressure  ratio,  which  produced  trailing  edge 
shocks  that  extended  to  the  exit  plane.  This  newer  version  allowed  for  a  supersonic  exit 
boxmdary  condition  where  the  exit  conditions  were  based  on  Giles’  characteristic  boundary 
conditions  [Ref.  11].  The  Baldwin-Lomax  turbulence  model  was  used  for  this  test  case  since 
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Figure  10.  Three-Dimensional  Multiple  Grid  (150  x  31  x  65) 
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a  'bug'  was  discovered  that  would  not  allow  the  code  to  properly  recognize  the  Cebeci-Smith 
turbulence  model  [Ref.  12].  The  code’s  author  has  since  corrected  the  program. 

Appendix  C  contains  an  example  namelist  input  file  used  to  obtain  a  flow  solution 
and  a  description  of  the  steps  required  to  run  the  code  on  the  Naval  Postgraduate  School 
(NPS)  and  National  Aerodynamic  Simulation  Program  (NAS)  Cray  supercomputers.  Output 
solution  file  (fort.3)  information  was  visually  examined  with  FAST  and  PLOT3D  graphics 
software  [Refs.  13  and  14].  Solution  residual  files  (fort.4),  blade  surface  pressure  files 
(fort. 7),  and  Mach  number  and  flow  angle  files  (fort.  7-11)  were  calculated  from  the  solution 
file  with  the  FORTRAN  programs  "pxy.f  and  "plane.f  (Appendix  G)  for  final  graphical 
representation. 
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IV.  RESULTS  AND  DISCUSSION 


A.  BLADE  MIDSPAN  SURFACE  PRESSURE  MEASUREMENTS 

Experimental  blade-surface  pressure  measurements  were  averaged  and  a  maximiun 
repeatability  error  was  calculated.  Each  blade  static  port  reading  was  non-dimensionalized  by 
the  upstream  stagnation  pressure  (Pq).  In  order  to  provide  an  experimental  comparison  with 
Reference  6  for  the  0.6815  pressure  ratio,  experimental  repeatability  was  confirmed  by 
conducting  six  runs;  three  runs  each  on  separate  days  (Tables  D1  and  D2  of  Appendix  D). 
The  maximum  repeatability  error  was  0.87%.  All  other  mean  pressure  data  were  based  on  a 
three  run  average  at  each  pressure  ratio  and  resulted  in  a  maximum  repeatability  error 
ranging  from  0.66%  at  the  lowest  pressure  ratio  down  to  0.22%  at  the  highest  pressure  ratio. 
Table  1  contains  each  midspan  static  port  location  and  non-dimensional  pressure 
measurement  for  all  experimental  pressure  ratios.  Tables  D3  through  D6  of  Appendix  D 
contain  pressure  data  for  pressure  ratios  of  0.5070,  0.6041,  0.8077,  and  0.9054  respectively. 
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0.9878 
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10 
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0.9600 
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11 
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0.8853 

0.8874 

0.9015 
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Table  1.  Non-Dimensional  Midspan  Surface  Pressure  (P/Po) 


B.  LASER-DOPPLER  VELOCIMETRY  MEASUREMENTS 

Two-dimensional  LDV  measurements  were  obtained  to  an  approximate  depth  of  3.34 
mm  through  a  1.0922  mm  (0.043  in)  diameter  hole  located  one-half  axial  chord 
downstream.  All  LDV  data  were  acquired  at  a  pressure  ratio  of  0.9054  and  are  tabulated  in 
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Tables  El  through  E5  of  Appendix  E.  AU  data  resulted  from  programmed  editing  of 
histograms  by  the  FIND  software  and  data  outside  two  standard  deviations  were  discarded. 

The  axial  velocity  (Vz)  and  tangential  velocity  (Vtheta)  were  combined  to  form  the 
total  velocity  (Vtotal)  shown  in  Figure  11.  Flow  velocities  are  seen  to  decrease  toward  the 
case  wall.  The  boundary  layer  was  distorted  due  to  secondary  flows  and  wake  and  comer 
vortices  that  forni  wdthin  the  blade  passage.  Two-dimensional  endwall  flow  angle  was  defined 
as  the  arc  tangent  of  the  tangential  velocity  divided  by  the  axial  velocity  and  is  shown 
graphically  in  Figure  12.  Periodicity  is  evident  over  11.6  degrees  (31  blades). 

Turbulence  intensities  were  calculated  with  respect  to  the  maximum  downstream  exit 
velocity  (Vexit)-  Figures  13  and  14  show  turbulence  intensity  in  both  the  tangential  and  axial 
directions.  The  tangential  turbulence  intensity  is  seen  to  be  higher,  possibly  due  to  the  steep 
gradient  of  the  Vtheta  mean-flow  profile  in  the  radial  direction. 

The  0.9054  pressure  ratio  provided  excellent  seeding  conditions  and  ATC  vibrations 
were  low.  Data  rates  ranged  from  approximately  300  samples  per  second  at  the  3.34  mm 
depth  to  150  samples  per  second  at  the  0.01  mm  depth.  Seeding  material  slowly  accumulated 
inside  the  lower  portion  of  the  optical  access  hole  and  interfered  with  LDV  data  acquisition. 
Occasionally,  the  atomizer  was  secured  and  seed  material  allowed  to  disperse.  Hole  alignment 
and  laser  power  (1.5  Watts)  were  periodically  verified  as  these  could  drift  due  to  temperature 
changes.  Endwall  measurement  techniques  did  not  provide  pressure  equalization  across  the 
optical  access  hole.  During  one  exploratory  mn  a  minimum  pressure  ratio  of  0.80  was 
achieved  before  the  data  rates  deteriorated  to  unacceptable  levels. 

A  random  mode  comparison  resulted  in  an  average  repeatability  difference  of  2.9%, 
0.6%,  6.4%,  and  7.0%  for  velocity,  flow  angle,  tangential  turbulence  intensity,  and  axial 
turbulence  intensity  respectively.  Figures  15  through  23  graphically  depict  velocity,  flow 
angle,  and  turbulence  intensity  random  mode  repeatability  data  for  wake  positions  of  +7,  0, 
and  -8  degrees.  A  random  and  coincidence  mode  comparison  resulted  in  an  average 
maximum  difference  of  0.4%,  0.2%,  4.1%,  and  6.3%  for  velocity,  flow  angle,  tangential 
turbulence  intensity,  and  axial  turbulence  intensity  respectively.  Tables  FI  and  F2  of 
Appendix  F  contain  LDV  repeatability  data. 
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Figure  12.  Two  Component  LDV  End  wall  Flow  Angle 
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Turbulence  Intensity  (%)  3  Turbulence  Intensity  (%) 


13.  Two  Component  LDV  EndwaU  Turbulence,  Vtheta 


Figure  14.  Two  Component  LDV  Endwall  Turbulence,  Vz 


Figure  16.  LDV  End  wall  Flow  Velocity  Repeatability  At  0  Degrees  Wake  Position 


LDV  Endwall  Flow  Velocity 
Repeatability  (-8°  Wake) 


Figure  17.  LDV  Endwall  Flow  Velocity  Repeatability  At  -8  Degrees  Wake  Position 
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Figure  18.  LDV  Endwall  Flow  Angle  Repeatability  At  +7  Degrees  Wake  Position 


Figure  19.  LDV  Endwall  Flow  Angle  Repeatability  At  0  Degrees  Wake  Position 


Figure  20.  LDV  Endwall  Flow  Angle  Repeatability  At  -8  Degrees  Wake  Position 
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Figure  22  LDV  Endwall  Turbulence  Repeatability  At  0  Degrees  Wake  Position 


LDV  Endwall  Turbulence  Intensity 
Repeatability  (-8°  Wake) 
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Figure  23.  LDV  Endwall  Turbulence  Repeatability  At  -8  Degrees  Wake  Position 
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C.  NUMERICAL  COMPARISON 

Computational  Fluid  Dynamic  (CFD)  solutions  were  run  with  the  same  pressure  ratios 
used  in  the  experiments.  Appendix  G  briefly  describes  the  computational  data  reduction 
process  and  programs  utilized.  Graphical  output  of  the  numerical  results  was  useful  in 
obtaining  a  qualitative  understanding  of  the  flow  characteristics. 

1.  Blade  Midspan  Surface  Pressures 

Reference  6  predicted  hub,  midspan  and  tip  surface  pressures  at  a  pressure  ratio  of 
0.68.  Figure  24  shows  a  comparison  between  numerical  and  experimental  blade  surface 
pressures  at  midspan  for  a  pressure  ratio  of  0.6815.  The  suction  surface  curve  suggests  that 
the  position  of  the  blade  passage  throat  was  at  0.8  of  axial  chord  and,  in  fact,  the  throat  was 
located  between  static  ports  five  and  six.  The  minimum  and  maximum  values  at  static  ports 
five,  six,  and  seven  may  be  due  to  the  boundary  layer  interaction  at  the  blade’s  blunt  trailing 
edge.  Figures  25  through  28  show  a  comparison  between  numerical  and  experimental  blade 
midspan  surface  pressures  at  pressure  ratios  of  0.5070,  0.6041,  0.8077  and  0.9054 
respectively.  Comparison  of  the  blade  surface  measurements  and  the  computational  results 
generally  show  excellent  agreement. 


0.6815  Pressure  Ratio 


Figure  24.  P/Pq  vs.  x/c  for  0.6815  Pressure  Ratio 
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P/PO  §  P/PO 


0.5070  Pressure  Ratio 


i  25.  P/Po  vs.  x/c  for  0.5070  Pressure  Ratio 


0.6041  Pressure  Ratio 


Figure  26  P/Pq  vs.  x/c  for  0.6041  Pressure  Ratio 
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27.  P/Po  vs.  x/c  for  0.8077  Pressure  Ratio 


0.9054  Pressure  Ra. 
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Figure  28.  P/Pq  vs.  x/c  for  0.9054  Pressure  Ratio 


2.  Mach  Number  and  Flow  Angle 

Two-dimensional  Mach  number  and  flow  angle  comparisons  were  made  between 
numerical  and  experimental  results  for  the  surveys  at  each  of  the  seven  radial  positions. 
Figures  29  through  32  compare  Mach  number,  while  Figures  33  through  36  compare  flow 
angle  at  four  selected  radial  positions.  The  remaining  comparisons  are  included  in  Appendix 
H.  All  comparisons  were  based  on  a  circumferential  match  between  the  experiment  and 
computation,  by  shifting  the  latter  to  coincide  with  the  measured  profile.  The  specific 
matching  was  based  on  the  maximum  Mach  number  at  the  deepest  radial  position  (3.34  mm). 
This  matching  was  kept  constant  for  the  circumferential  comparisons  at  all  other  spanwise 
locations.  To  depict  flow  periodicity  and  equivalently  compare  with  experimental  wake 
positions,  the  numerical  solutions  were  repeated  over  one  and  one-half  blade  passages. 

Computed  Mach  numbers  differed  from  the  experimental  values  by  an  average 
12.7%  between  the  surveys  at  0.18  mm  and  3.34  mm.  The  difference  increased  to  70.8% 
near  the  endwaU,  possibly  due  to  a  combination  of  LDV  data  velocity  biasing  and  insufficient 
grid  resolution.  Velocity  biasing  was  estimated  to  cause  an  approximate  11.0%  velocity 
increase  in  the  endwall  region.  The  numerical  solution  at  the  endwaU  (0.01  mm)  was  based 
on  the  last  k  grid  point  in  the  radial  direction.  Computed  flow  angles  surprisingly  only  varied 
from  the  experimental  values  by  an  average  3.2%  (2.3  degrees).  This  was  in  contrast  to  the 
measurements  in  Reference  6  (Prat  =  0.68)  which  showed  poor  comparison  between 
measured  and  computed  flow  angles  across  the  wake  at  90  percent  span. 

The  question  of  radial  spatial  error  with  the  probe  volume  needed  to  be  determined. 
The  probe  volume  length  and  diameter  (Figure  37)  was  calculated  as  1.56  mm  and  0.11  mm 
respectively.  The  probe  volume  had  approximately  30  fiinges  across  its  minor  axis  and  after 
data  processing  an  average  of  13  fringe  crossings  constituted  a  valid  Doppler  burst.  Since  the 
IFA-750  digital  burst  correlator  automatically  centered  the  valid  signal  within  the  probe 
volume,  an  effective  probe  volume  with  length  1.39  mm  and  diameter  0.04  mm  was  obtained. 
Figure  38  depicts,  to  scale,  the  radial  survey  resolution  that  results  from  endwall  depth  and 
effective  probe  volume  dimensions.  The  decrease  in  relative  spatial  resolution  close  to  the 
endwall  also  contributed  to  the  lack  of  agreement  between  the  Mach  number  profiles  (Figure 
32). 
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29.  Mach  Number  Comparison  At  3.34  mm  Depth 
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Figure  30.  Mach  Number  Comparison  At  0.89  mm  Depth 
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j  31.  Mach  Number  Comparison  At  0.18  mm  Depth 
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Figure  32.  Mach  Number  Comparison  At  0.01  mm  Depth 
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33.  Wake  Flow  Angle  Comparison  At  3.34  mm  Depth 


Figure  34.  Wake  Flow  Angle  Comparison  At  0.89  mm  Depth 
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35.  Wake  Row  Angle  Comparison  At  0.18  mm  Depth 
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Figure  38.  Endwall  Measurement  Schematic 
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D.  NUMERICAL  RESULTS 


1.  Convergence  History 

Solutions  at  pressure  ratios  of  0.5070,  0.6041  and  0.6815  fully  converged  after 
approximately  1,200  iterations.  Figure  39  shows  density  residuals  at  0.6815  pressure  ratio 
after  10,200  iterations.  Solutions  for  pressure  ratios  of  0.8077  and  0.9054  were  initially 
obtained  after  1,200  iterations,  but  convergence  did  not  occur  (i.e.  three  orders  of  magnitude 
reduction  of  the  residuals)  until  approximately  3,000  iterations.  Figures  II  through  15  of 
Appendix  I  contain  density-residual  convergence  history  for  all  examined  pressure  ratios. 


Figure  39.  0.6815  Pressure  Ratio  Convergence  History  (10,200  Iterations) 


2.  Turbulence  Models 

At  a  pressure  ratio  of  0.8077,  the  code  was  run  with  the  following  turbulence  models; 
RVC3D  version  206  Cebeci-Smith 

RVC3D  version  206  Baldwin-Lomax 

RVC3D  version  208  Baldwin-Lomax 

so  as  to  assess  the  effect  of  different  turbulence  models  on  the  solution.  In  aU  three  cases,  the 
code  produced  similar  solutions. 
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3.  Mach  Number  Contours 

Figures  40  and  41  display  exit  plane  Mach  number  contours  at  pressure  ratios  of 
0.9054  and  0.5070  respectively.  The  Mach  number  contours  between  25%  and  75%  span 
provide  an  indication  of  where  the  flow  characteristics  were  well  established.  Boundary  layers 
appear,  as  expected,  to  be  thin  near  the  endwall  region  and  thicker  near  the  hub.  Spanwise 
wake  curvamre,  due  to  secondary  flows,  is  also  apparent  in  the  two  figures  with  Figure  40 
showing  the  greater  curvature. 

Figure  42  displays  midspan  Mach  number  contours  for  the  blade  passage  at  0.5070 
pressure  ratio.  The  flow  was  strongly  accelerated  along  the  suction  side.  The  acceleration  of 
the  flow  to  high  local  Mach  numbers  resulted  in  shocks  at  the  throat  and  trailing  edge.  The 
wake  was  a  merging  of  the  pressure  and  suction  side  boundary  layers.  The  trailing-edge 
shock  extended  across  the  wake  from  the  adjacent  (upper)  blade,  out  to  the  exit  plane  of  the 
grid.  The  shock  interaction  with  the  wake  caused  the  wake  to  narrow  down,  resulting  in  a 
wake  which  was  repeatedly  diffused  and  then  coalesced.  The  shock  strength  decreased  away 
from  the  trailing  edge  as  a  result  of  this  interaction  with  the  wake,  which  could  be  both 
physical  and  computational  in  nature.  Because  of  the  relatively  good  agreement  between 
experiment  and  computation  of  the  blade  surface  pressure  field,  it  was  felt  that  the  solution 
was  realistic.  Appendix  J  includes  program  listings  for  a  graphical  display  of  multiple  grids 
and  corresponding  numerical  solutions  for  use  with  PLOT3D. 
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Figure  40.  Exit  Plane  Mach  Contours  (0,9054  Pressure  Ratio) 
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Figure  42.  Midspan  Mach  Contours  (0.5070  Pressure  Ratio) 


V.  CONCLUSIONS  AND  RECOMMENDATIONS 


An  experimental  and  computational  investigation  was  conducted  of  the  flow  through 
an  annular  turbine  cascade.  Computational  pressure  ratios  (Prat)>  defined  as  the  downstream 
hub-static  pressure  (Phub)  divided  by  the  upstream  stagnation  pressure  (Pq),  were  chosen  to 
coincide  with  experimental  test  conditions.  Downstream  computational  and  experimental 
measurement  locations  were  chosen  to  coincide  to  provide  a  proper  comparison. 

An  existing  annular  turbine  cascade  was  modified  to  provide  a  capability  to  measure 
blade  midspan  surface  pressures.  Experimental  measurements  of  blade  midspan  surface 
pressures  were  obtained  and  favorably  compared  with  numerical  predictions  over  a  range  of 
subsonic  and  transonic  conditions.  Measurements  were  repeatable  with  a  maximum 
repeatability  error  of  less  than  0.87%.  Numerical  solutions  converged  on  pressure  within 
approximately  1,200  iterations. 

The  minimum  and  maximum  pressure  values  on  the  suction  surface  close  to  the 
trailing  edge  should  be  investigated  further  to  confirm  the  influence  of  the  relatively  blunt 
trailing  edge.  An  inlet  boundary  layer  survey,  requiring  a  rig  modification,  should  be 
performed  to  supplement  data  in  Reference  6  and  refine  the  inlet  boundary  conditions  to  the 
computational  model. 

To  supplement  blade  midspan  surface  pressures,  blade  pressure  taps  could  be  installed 
at  other  spanwise  locations  to  include  the  hub,  tip,  and  possibly  endwaU  locations.  A  tip  gap 
could  be  created  and  a  computational  tip  gap  model  introduced  to  investigate  the  endwaU 
region.  Blade  tips,  and  the  adjacent  endwaU,  could  be  instrumented  for  pressure  readings  and 
a  numerical  and  experimental  comparison  conducted  to  obtain  insight  into  tip  leakage  flows. 
In  addition,  an  LDV  optical-access  window  could  be  instaUed  over  the  tip  region  to  measure 
tip-gap  flow  characteristics. 

Successful  measurements  of  two-dimensional  velocity,  flow  angle,  and  turbulence 
intensity  through  a  1.0922  mm  diameter  casing-access  hole  were  demonstrated.  All  LDV 
data  were  acquired  one-half  axial  chord  downstream  over  seven  radial  positions  and  17 
separate  peripheral  displacements.  Experimental  flow  velocities  were  increased  from  those  of 
Reference  6  mainly  due  to  the  availability  and  use  of  a  digital  burst  correlator  in  place  of  a 
signal  processor.  Although  all  LDV  data  reported  herein  were  acquired  at  0.9054  pressure 
ratio,  measurements  were  attempted  at  0.80  pressure  ratio  and  some  data  were  obtained; 
however,  low  data  rate  and  endwaU  hole  vibrations  prevented  a  complete  survey. 
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The  comparison  of  computational  and  experimental  Mach  number  and  flow  angle 
yielded  an  average  12.7%  difference  in  Mach  number  and  3.2%  (2.3  degrees)  difference  in 
flow  angle.  AH  experimental  measurements  were  repeatable  to  an  average  uncertainty 
ranging  from  0.2%  to  7.0%  for  velocity,  flow  angle,  and  turbulence  intensity  in  both  random 
and  coincidence  data  processor  modes.  All  numerical  and  experimental  comparisons  were 
based  on  a  circumferential  match  of  maximum  Mach  number  at  the  80.9%  (3.34  mm)  span 
location.  An  investigation  of  other  circumferential  matching  options  should  be  conducted  to 
explore  the  effect  on  the  degree  of  agreement  (i.  e.  minimum  Mach  number  circumferential 
matching). 

LDV  alignment  and  radial  survey  procedures  were  devised  and  repeated.  The  probe 
volume  passed  cleanly  through  the  optical  access  hole  and  radial  surveys  were  conducted  to 
provide  complete  coverage  for  data  recording.  Probe  volume  dimensions  were  further 
refined  to  aid  in  estimates  of  radial-position  accuracy.  Introduction  of  a  second  fiber-optics 
probe  for  three  component  measurements  would  decrease  the  effective  probe  volume  and 
increase  radial  position  accuracy.  The  ability  to  resolve  endwaU  flow  characteristics  and 
eventually  flow  tip  losses  appears  to  be  promising. 

The  present  hardware  did  not  provide  pressure  equalization  across  the  optical  access 
hole.  A  rig  modification  is  required  to  equalize  the  pressure  at  the  measurement  location  and 
eliminate  any  influence  the  opening  might  have  on  the  flow  characteristics.  Also,  the  coarse 
peripheral  (wake)  positioning  mechanism  needs  modification  to  allow  more  precise 
adjustment.  Experimental  peripheral  positioning  uncertainties  may  contribute  to  differences 
in  the  compared  data. 

Future  efforts  with  computational  fluid  dynamics  should  include  modeling  the  inlet 
bellmouth  and  center  body,  and  conducting  a  turbulence  model  sensitivity  analysis  in  the 
endwall  region.  Increased  grid  resolution  in  the  endwaU  region  would  be  beneficial.  The 
experimental  99.9%  span  location  coincided  with  the  final  k  grid  point  and  may  account  for 
the  70.8%  difference  between  the  experimental  and  computed  Mach  numbers.  Finally,  a 
detailed  investigation  into  transonic  and  sonic  flow  conditions  could  be  conducted  to  examine 
the  blade-passage  and  trailing-edge  shock  locations  and  subsequent  interactions,  both 
computationally  and  experimentally. 
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APPENDIX  A.  PRESSURE  DATA  ACQUISITION 


Pressure  data  acquisition  was  accomplished  using  a  Hewlett-Packard  9000  computer 
system  with  the  following  program: 
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iHLENAME:  BILLSP 
IDISK  LABEL:  ’"/AUSTIN" 

ILAST  MODIFIED  10/27/94  BY  LT  BILL  DONOVAN 

ITHIS  PROGRAM  RECORDS  AND  REDUCES  MEASURED  PRESSURES  FROM  A 
ISCAMVALVE  CONNECTED  TO  11  STATIC  PRESSURE  PORTS  ON  ATC 
!  VARIABLES  /  PORT  CONNECTIONS... 

!  V=DESIRED  S/V 

!  A=PRESENT  S/V  PORT 

!  C=CHANNEL 

!  Voltage  (  )=DVM  VOLTAGE  READING  FOR  EACH  PORT 
!  S=SCANNER  NUMBER 

!  Pinhg=BAROMETRIC  PRESSURE  IN  INHG 

!  Pamb=BAROMETRIC  PRESSURE  (psi) 

!  AMBIENT  PRESSURE...#! 

!  CALIBRATION  PRESSURE...#2 

!  REFERENCE  UPSTREAM  TOTAL  PRESSURE... PORT  #4 

!  INNER  HUB  STATIC  PRESSURE.. .PORT  #5 

!  UPSTREAM  STATIC  PRESSURE...PORT  #6 

!  ATC  PASSAGE  PORTS  #1-#1 1...S/V  PORTS  #9-#19 

!  Prat=  INNER  HUB  STATIC/  REFERENCE  UPSTREAM  TOTAL 
!  Tstag=UPSTREAM  STAGNATION  TEMPERATURE  IN  DEGREES  F 

OPTION  BASE  1 
CLEAR  SCREEN 
PRINTER  IS  CRT 

DISP  "PLEASE  WAIT  WHILE  RESETTING  SCANIVALVE  TO  PORT  #1" 

PRINT 

V=5 

A=1 

GOSUB  Read 
DIM  Y(20) 

DIM  Voltage(50) 

DIM  P(50) 

INPUT  "ENTER  MONTH,  DAY,  YEAR  (I.E.02, 04,  94)",  Y(3),  Y(5),  Y(7) 

INPUT  "ENTER  RUN  #: ",  Y(9) 

PRINTER  IS  711 

PRINT  USING  "K,  DD,  "  "  /  "  ",  DD, "  "  /  "  ",  DD",  "DATE  OF  RUN:  ",  Y(3),  Y(5),  Y(7) 
PRINT  "DATA  RUN  ",  Y(9) 

PRINT  USING  "A  5X,  K,  DDD,  /,  2X,  K,  15X,  K,  17X,  K”,  "SCANIVALVE  #  ",  V, 
"Port",  "Volts",  "Psi" 

PRINTER  IS  CRT 

PRINT  "  ZERO  DVM  ON  SCANIVALVE  #5,  PORT  #1 ..." 

PRINT 

PRINT  "  SET  20  inhg  FOR  SCANIVALVE  CALIBRATION..." 

PRINT 

PRINT "  ADJUST  SPAN  FOR  DVM  =  0.009823  VOLTS  ON  PORT  #2..." 
PRINT 

PRINT  "ENSURE  DESIRED  PRESSURE  RATIO  IS  SET..." 

INPUT  "ENTER  BAROMETRIC  PRESSURE  IN  INHG",  Pinhg 
INPUT  "ENTER  HIGHEST  SCANIVALVE  PORT  #  NEEDED  ”,  Nports 
INPUT  "ENTER  UPSTREAM  STAGNATION  TEMP.  IN  DEG.  F  ",  Tstag 
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Pamb=Pinhg*  (0.491 1541) 

PRINT  "****  PRESS  CONTINUE  WHEN  READY  TO  TAKE  DATA***” 
PAUSE 

I****************************************************** 

iRECORD  THE  PRESSURE  DATA  FOR  THIS  RUN 

OUTPUT  722;  "FIRITIMOZI"  !  SETS  UP  HP3456A  DVM 
V=5  ISCANIVALVE  #5  DESIGNATED 

S=1  IHP3495A  SCAN  #1  DESIGNATED 

FOR  A=1  TO  Npoits 
GOSUB  Read 
WAIT  2.0 

OUTPUT  701;  "C"  ICLEARS  SCANNER  #1 

OUTPUT  701  USING  "DDD";  V+9  !SCAN  SET  TO  READ  SA^ 
TRIGGER  722 
ENTER  722;  Voltage(A) 

P(A)=(Voltage(A)  -  VoItage(I))*l(X)0+Pamb 
PRINTER  IS  CRT 

PRINT  USING  "2X,DDD,7X,7DD.DDDDDDD,8X,7DD.DDDD";  A,  VoItage(A),  P(A) 
PRINTER  IS  711 

PRINT  USING  "2X,DDD,7X,7DD.DDDDDDD,8X,7DD.DDDD";  A,  Voltage(A),  P(A) 
NEXT  A 
Prat=P(5)/P(4) 

PRINT  "Prat  =",  Prat 

PRINT  "UPSTREAM  STAGNATION  TEMPERATURE  WAS  (DEG  F)",  Tstag 
PRINT  "INPUT  BAROMETRIC  PRESSURE  IN  (inhg)  WAS  ",  Pinhg 
j*****=ic****suBROUTINE  TO  POSITION  AND  READ  SfV  PORTS********* 
Read:  OUTPUT  707  USING  "#,K";V 
PO=SPOOL(707) 

L=BINAND(P0,15) 

T1=ROTATE(PO,4) 

M1=BINAND(T1,7) 

P(A)=10*M1+L 
CLEAR  707 

IF  P(A)=A  THEN  Finish 
OUTPUT  701;  "C" 

OUTPUT  701  USING  "DDD";V-I 
OUTPUT  701;  "C" 

WAIT  0.1 
GOTO  Read 
Finish:  RETURN 
END 
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The  following  table  relates  each  Scanivalve  port  to  its  respective  pressure 
measurement. 


Scanivatv«  Port  Number 

Pressiffid'Measured 

1 

Ambient  Pressure 

2 

Calibration  Pressure  (set  at  20  inches  hg) 

3 

Not  Used 

4 

Upstream  Total  Pressure  (Po) 

5 

Downstream  Hub  Static  Pressure  (Phub) 

6 

Upstream  Static  Pressure  (Pstatic) 

7 

Not  Used 

8 

Not  Used 

9 

Blade  Static  Port  #1  (at  leading  edge) 

10 

Blade  Static  Port  #2  (suction  side) 

11 

Blade  Static  Port  #3  (suction  side) 

12 

Blade  Static  Port  #4  (suction  side) 

13 

Blade  Static  Port  #5  (suction  side) 

14 

Blade  Static  Port  #6  (suction  side) 

15 

Blade  Static  Port  #7  (suction  side) 

16 

Blade  Static  Port  #8  (at  leading  edge) 

17 

Blade  Static  Port  #9  (pressure  side) 

18 

Blade  Static  Port  #10  (pressure  side) 

19 

Blade  Static  Port  #1 1  (pressure  side) 

Table  Al.  Pressure  Data  Acquisition  Connections 
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APPENDIX  B.  GRID  GENERATION  INPUT  FILE 


Grid  generation  was  completed  on  Silicon  Graphics  Indigo  II  workstations.  The 
following  TCGRID  namelist  file  was  utilized  to  obtain  the  final  grid. 


&naml  im=150  jm=31  km=65  itl=20  icap=18  k2d=3  merid=0  &end 
&nam2  nle=16  nte=10  dsle=.018  dste=.(X)3  dshub=.(X)04  dstip=.0004 
dswte=.001  dswex=.060  dsthr=l.  dsmin=.0004  dsmax=.025 
dsra=.45  rcom=.098  &end 

&nam3  iterm=100  idbg=0  0  0  0  0  0  0  0  aabb=.5  ccdd=.45  &end 
&nam4  zbc=-1.5  -1.5  0.5  -1.5  -1.5  .5 

rbc=  3.895  3.895  3.895  4.585  4.585  4.585  &end 
TRANSONIC  TURBINE’ 

2  2 


-.975 

0.0 

3.895 

3.895 

-.975 

0.0 

4.585 

4.585 

2 

51 

-7.8999996E-03  -9.9200001E-03  -1.2000000E-02  -1.4080000E-02 
- 1.61 0000 lE-02  -1.7999999E-02  -1.9710001E-02  -2.1190001E-02 
-2.2390001E-02  -2.3280000E-02  -8.1040002E-02  -0.1503800 
-0.2308900  -0.3219400  -0.4228500  -0.5328500  -0.6511000 
-0.7893000  -0.9275000  -0.9410600  -0.9527700  -0.9622700 
-0.9692700  -0.9735600  -0.9750000  -0.9735600  -0.9692700 
-0.9622700  -0.9527700  -0.9410700  -0.9275000  -0.9124900 
-0.8250000  -0.7200000  -0.6000000  -0.4560000  -0.3440000 
-0.2730000  -0.2200000  -0.1820000  -0.1160000  -5.9000000E-02 
-1.7000001E-02  O.OOOOOOOE+00  -1.8220000E-04  -7.2359998E-04 
-1.6100000E-03  -2.8100000E-03  -4.2900001E-03  -6.0000001E-03 
-7.8999996E-03  -0.2349885  -0.2351271  -0.2351733  -0.2351271 
-0.2349885  -0.2347600  -0.2344519  -0.2340719  -0.2336329 

-0.2331451  -0.1980565  -0.1643980  -0.1324210  -0.1023671 

-7.4464694E-02  -4.8926830E-02  -2.5946086E-02  -4.6726577E-03 
1.6618744E-02  1.9057767E-02  2.2061616E-02  2.5545571E-02 
2.9399229E-02  3.3504494E-02  3.7740692E-02  4.1976891E-02 
4.6082158E-02  4.9935814E-02  5.3417202E-02  5.6423619E-02 
5.8862645E-02  6.0659818E-02  6.9319643E-02  7.1887039E-02 
6.7522466E-02  5.1347882E-02  2.5673941E-02  O.OOOOOOOE+00 
-2.5673941E-02  -5.1347882E-02  -0.1026958  -0.1540437 
-0.2053915  -0.2320924  -0.2326265  -0.2331451  -0.2336329 

-0.2340719  -0.2344519  -0.2347600  -0.2349885  3.895000 

3.895000  3.895000  3.895000  3.895000  3.895000  3.895000 

3.895000  3.895000  3.895000  3.895000  3.895000  3.895000 

3.895000  3.895000  3.895000  3.895000  3.895000  3.895000 

3.895000  3.895000  3.895000  3.895000  3.895000  3.895000 

3.895000  3.895000  3.895000  3.895000  3.895000  3.895000 

3.895000  3.895000  3.895000  3.895000  3.895000  3.895000 

3.895000  3.895000  3.895000  3.895000  3.895000  3.895000 

3.895000  3.895000  3.895000  3.895000  3.895000  3.895000 

3.895000  3.895000  -7.8999996E-03  -9.920000 lE-03 
-1.2000000E-02  -1.4080000E-02  -1.6100001E-02  -1.7999999E-02 
-1.9710001E-02  -2.1190001E-02  -2.2390001E-02  -2.3280000E-02 
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-8.1040002E-02  -0.1503800  -0.2308900  -0.3219400  -0.4228500 
-0.5328500  -0.6511000  -0.7893000  -0.9275000  -0.9410600 
-0.9527700  -0.9622700  -0.9692700  -0.9735600  -0.9750000 
-0.9735600  -0.9692700  -0.9622700  -0.9527700  -0.9410700 
-0.9275000  -0.9124900  -0.8250000  -0.7200000  -0.6000000 
-0.4560000  -0.3440000  -0.2730000  -0.2200000  -0.1820000 
-0.1160000  -5.9000000E-02  -1.7000001E-02  O.OOOOOOOE-i-00 
-1.8220000E-04  -7.2359998E-04  -1.6100000E-03  -2.8100000E-03 
-4.2900001E-03  -6.000000 lE-03  -7.8999996E-03  -0.1996249 
-0.1997426  -0.1997819  -0.1997426  -0.1996249  -0.1994308 
-0.1991690  -0.1988462  -0.1984733  -0.1980589  -0.1682508 
-0.1396576  -0.1124929  -8.6961828E-02  -6.3258447E-02 
-4.1563794E-02  -2.204 1440E-02  -3.9694658E-03  1.4117776E-02 
1.6189748E-02  1.8741548E-02  2.1701200E-02  2.4974918E-02 
2.8462378E-02  3.2061070E-02  3.5659760E-02  3.9147221E-02 
4.2420939E-02  4.5378406E-02  4.7932386E-02  5.0004359E-02 
5.1531076E-02  5.8887679E-02  6.1068702E-02  5.7360962E-02 
4.3620501E-02  2.1810250E-02  O.OOOOOOOE+00  -2.1810250E-02 
-4.3620501E-02  -8.7241001E-02  -0.1308615  -0.1744820 
-0.1971647  -0.1976183  -0.1980589  -0.1984733  -0.1988462 
-0.1991690  -0.1994308  -0.1996249  4.585000  4.585000  4.585000 

4.585000  4.585000  4.585000  4.585000  4.585000  4.585000 

4.585000  4.585000  4.585000  4.585000  4.585000  4.585000 

4.585000  4.585000  4.585000  4.585000  4.585000  4.585000 

4.585000  4.585000  4.585000  4.585000  4.585000  4.585000 

4.585000  4.585000  4.585000  4.585000  4.585000  4.585000 

4.585000  4.585000  4.585000  4.585000  4.585000  4.585000 

4.585000  4.585000  4.585000  4.585000  4.585000  4.585000 

4.585000  4.585000  4.585000  4.585000  4.585000  4.585000 


Figure  Bl.  Blade  Geometry  From  Ref.  [6] 


46 


APPENDIX  C.  RVC3D  SAMPLE  INPUT  AND  EXECUTION 


The  following  input  namelist  example  was  for  a  0.6815  pressure  ratio  with  subsonic 
exit  boundary  conditions  and  Cebeci-Smith  turbulence  modeling. 


TRANSONIC  TURBINE  Annular  Cascade' 

&nll  im=150  jm=31  km=65  itl=20  iil=67  «feend 

&nl2  cfl=5.5  aviscl=0.0  avisc2=0.0  avisc4=0.50  ivdt=l  nstg=4  itmax=1200 
irs=l  epi=0.50  epj=0.60  epk=0.60  &end 
&nl3  ibcin=3  ibcex=3  isymt=0  ires=10  icmt=50 
iresti=0  iresto=l  ibcpw=0  iqin=0  &end 
&nl4  emxx=0.13  emty=0.0  emrz=0.0  expt=0.0  prat=0.6815  ga=1.4 
om=0.000000  igeom=l  alex=-67.0  &end 
&nl5  ilt=3  tw=1.00  renr=6.651e6  pmr=.7  prtr=.9  vispwr^.666666 

srtip=0.0  cmutm=14.  jedge=31  kedge=20  iltin=2  dblh=0.0048  dblt=0.0178  &end 
&nl6  iol=l  io2=165  oar=0.  ixjb=0  njo=l  nko=3 
jo=l  ko=5  11  16  &end 


Runs  on  the  NPS  Cray  computer  took  approximately  12  hours  for  a  1,200  iteration 
convergence  and  were  accomplished  with  the  command: 

>qsub  myjob.nqs 

where  the  file  myjob.nqs  was  as  follows: 
myjob.nqs 

#QSUB  -q  prem  -IT  50000  -IM  9Mw 
cd  /dl/whdonova 

./exec.rvc3d  <  rvc3d.in  >  rvc3d.out 


Runs  on  the  NASA  Ames  Research  Center  Cray  computer  took  approximately  1  hour 
for  a  1,200  iteration  convergence  and  were  accomplished  with  the  command: 

>qsub  -Im  lOMw  -It  14400  rvc3d.com 
where  the  file  rvc3d.com  was  as  follows: 
rvc3d.com 

cd  An/vn/v3/wdonovan 
./exec.rvc3d  <  rvc3d.in  >  rvc3d.out 
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APPENDIX  D.  BLADE  MIDSPAN  SURFACE  PRESSURE  DATA 


m?2ijgg2j| 

tagin; 


Po  (Ref.  Upstream  Stagnation) 


P  Static  Port  1  (@  leading  edge) 


ll£E&3£S£l 

m 

IXSESLSE&El 


rffljgjafi 


Table  Dl.  0.6815  Pressure  Ratio  Series  #1 
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Ityrnggiw—— 


Po  (Ref.  Upsiream  Stagnation) 


i!i213E53^SI^5I^^^E]^3ES^ES^S23] 


P  up  static  (Ref.  Upstream  Static) 


P  Static  Port  1 «®  leading  edge) 


latfgnffjgffTiM 


P  Stauc  Port : 


P  Static  Port  4 


ig«n:Tagfft»a 


P  Static  Port 


P  Static  Port  9 


P  Static  Port  10 


Table  D3.  0.5070  Pressure  Ratio  Series 
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Pressw«  (p«i> 


nLyalve  rsrt 


1 


cattivatverort  Yaltte 


P  ambient 


P  calibration 


emaMM 


l£ill!£iflm^SS5II^ 

\i 


hub  exit  (Static  @  0.5c  downstream) 


P  up  static  (Ref.  Upstream  Static) 


iganxCTggTwrgagPtnn 


P  Static  Port  3 


tatic  Port 


ig.-ffnagCT>i 


P  Stauc  Port 


taticport 


P  Sialic  Port  8  ((2)  leading  edge) 


tatic  Port  9 


rsR? 
lawmgagmiPi 


P  Static  Port  11 


Table  D5.  0.8077  Pressure  Ratio  Series 
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APPENDIX  E.  LDVDATA 


^ach 


Number 


turb.  int  I  turb.  int  |  f 


[)6  I  M  exit 


T  exit  CK) 


0.9056 

0379 

299.4 


m 


ake  Position 


[)epth 


Number 


turb.  int.  I  turb.  int.  I  P  static  (mb 


a  exit  (m/s) 
a  stag  {mis] 


Table  El.  LDV  Data  (Wake  Positions  of  -8,  -7,  -6,  and  -5  Degrees) 
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Table  E2.  LDV  Data  (Wake  Positions  of  -4,  -3,  -2,  and  -1  Degrees) 
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Table  E3.  LDV  Data  (Wake  Positions  of  0,  1, 2,  and  3  Degrees) 
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Table  E5.  LDV  Data  (Wake  Position  of  8  Degrees) 


APPENDIX  F.  LDV  REPEATABILITY  DATA 


inadence 


mcidenceO  1.78  104.0  24.7  106.9  0.304  76.6  4.33  |  3.33 


27.1  110.4  0.314  75.8  3.39 


incidence  0  I  3.34 


351.15 


V  exit  (m/s)  131.23 


Wake  Position  1  Depth  ]  V-theta 


V  total  j  Mach  J  Flow  Angle 


umber 


1  turb.  int. 

turb.  int. 

00 

4.83 

mgiii 

mygyu 

mam 

myni^ 

Table  FI.  LDV  Repeatability  Data  (Wake  Positions  of  -8  and  0  Degrees) 
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|jS3i3 

mym 

ygjjgjgg 

V-z 

imgmm 

mmi 

1  (degrees) 

IQgl 

KQgl 

^Qgjl 

msi 

Uirb.  int 

turb.  int 

mfims 

mmm 

mm 

g^mi 

mmii 

mym 

mimy 

jnnmm 

igmmi 

mmi 

mmm 

mg| 

mmi 

1  24.2 

mmm 

mmmi 

m^^ii 

imm 

mgmm 

Prat 

0.9055 

0.18 

92.1 

jmm 

mn^m 

0.271 

751 

mmi 

mmm 

M  exit 

0379 

jgmi 

gmm 

24,4 

mmi 

imm 

7^.4 

mmm 

mmm 

mmmm 

299.4 

Coincidence  7 

imi 

^mm 

mmii 

mmi 

0.288 

— 7525 — 

mmm 

mmm 

BBSfifiSI 

Coincidence  7 

mm 

mmn 

mmi 

mmm 

mmm 

76.4 

mmm 

mmm 

mg^gom 

mmm 

1 

mm 

mmi 

mmm 

mmm 

mmm 

16,1 

mmmi 

mmm 

mmsim 

mm 

Wake  Position 

Igjgl 

V- theta 

mym 

V  total 

Mach 

t low  Angle 

V -theta 

V-2 

Tstag(°D 

95.0 

(degrees) 

BSEBI 

(m/s) 

(m/s) 

(m/s) 

mumm 

(degrees) 

turb.  int 

P  Static  (inhg) 

30.00 

mm 

m^m 

mmm 

mmi 

mmm 

74.6 

^mm 

^^m 

mym 

mm 

mmi 

mmm 

mmm 

my^m 

7T7 

mymm 

mmm 

Prat 

mimm 

mm 

mmi 

mmi 

mmm 

mmi 

7525 

mmm 

mmm 

mmm 

mgm 

mm 

Kyi 

mmm 

mmm 

mmm 

7515 

^mm 

mmm 

299.4 

mm 

m^m 

mmi 

101.1 

mmH 

76.1 

^mm 

imm 

a  exit  (m/s) 

34634 

mmi 

mgm 

24.4 

163.9 

mmi 

75^4 

mmm 

4.66 

a  Stag  (m/s) 

351.79 

mm 

mmi 

35!3 

mmm 

mmm 

753 

mmm 

mmm 

mgmm 

mmmi 

Table  F2.  LDV  Repeatability  Data  (Wake  Position  of  7  Degrees) 
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APPENDIX  G.  COMPUTATIONAL  DATA  REDUCTION 


Residual  (fort.4)  and  P/Pq  (fort. 7)  files  were  produced  from  the  solution  (fort.3)  files 
using  the  FORTRAN  program  pxy.f. 


Program  pxy.f 

c  pxy.f  reads  rvc3d  files  &  writes  ascii  files  for  gnuplot 
c  unit  1  =  input  xyz  file 
c  unit  3  =  input  q  file 

c  unit  7  =  output  blade  pressures  on  5  k-planes 

c  unit  4  =  output  residual  history 

parameter(ni=150,nj=3 1  ,nk=65) 

real  x(ni,nj,nk),y(ni,nj,nk),z(ni,nj^) 
real  qq(5,ni,nj,nk),resd(5000,5) 
real  pk(5),xk(5) 
dimension  kk(5) 

c  k-values  are  hard-wired  below  (hub,  25%,  50%,  75%,  tip) 
data  kk/2,27,33,39,64/ 

^Hc**:je**4e***5|c:ie*:>:****He******5ic******3jcJ|e*:Jc*5je:|eJjc5{c**:<c*****3le;ie4t*:jc**:Jc:ie:|e** 

c  read  grid  coordinates 
read(l,*)im,jm,km 

read(  1  ,*)(((x(i  J,k),i=l  ,im),j= 1  ,jm),k=  1  ,km), 

1  (((y(i.j.k),i=l  ,im),j=  1  ,jm),k=l  ,km), 

2  (((z(ij,k),i=l,im),j=ljm),k=l,km) 

c  read  restart  file 

read(3,*)imax,jmax,kmax 

read(3  ,*)fsmach,alpha,re,time 
c 

icheck=iabs(im-imax)+iabs(jm-jmax)+iabs(km-kmax) 

if(icheck.ne.O)then 

write(6,610)im,jm,km,imax,jmaxjcmax 

stop 

endif 

c 

read(3,*)((((qq(l,i,j,k),i=l  ,im),j=l  ,jm),k=  l,km),l=l  ,5) 
c 

c  additional  residual  data 

read(3,*)itl,iil,phdeg,ga,om,nres,dum,dum,dum,dum 
read(3,*)((resd(nr,l),ni^  1  ,nres),l=  1 ,5) 

c  ps/pO  output  to  unit  3 

itr=im+l-itl 

ggm=ga*(ga-l.) 

j=l 
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c  nonnalize  x  by  chord 
do  7 1=1,5 

k=kka) 

xmin=x(im^,j,k) 
xmax=x(itl,  j,k) 
do  5  i=itl,itr 

xniin=min(xmin,x(i,j,k)) 

5  xmax=max(xmax,x(i,j,k)) 
chord=xmax-xinin 
do  7  i=itl,itr 

7  x(i,j,k)=(x(i,j,k)-xmin)/chord 

write(7,305)(kka)a=l,5) 
do  20  i=itl,itr 
do  10 1=1,5 

lf=lflf(1) 

pk(l)=ggm*(qq(5,i,j,k)-.5*(qq(2,i,j,k)**2+qq(3,i,j,k)**2 

l+qq(4,i,jdc)**2Vqq(l,i,j,k)) 

10  xkO)=x(i,j,k) 

20  write(7,300)i,(xka),pk(l),l=l,5) 
c  residual  history  output  to  unit  4 

^Hc****:ic***4«*  ************  **********Ji«  *****************  ********** 

write(4,310)  l,(resd(l,l),l=l,5) 
do  30  j=2,nres 

it=10*a-l) 

30  write(4,310)it,(resd(j,l),l=l,5) 

^****Jie*******************5i«****5fc******5ic********sJf*****  ********** 

300  format(i5,10fl0.5) 

305  formate  k=',5(17x,i3)) 

310  fonnat(i5,5(lx,el0.3)) 

610  formate  *****  warning  *****7, 

1  '  im,  jm,  km,  read  from  input' ,3i5,’  do  not  match'y, 

2  ’  im,  jm,  km,  read  from  restart  file’,3i5) 
stop 

end 
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The  plane.f  program  from  Reference  6  was  modified  to  compare  two-dimensional 
experimental  and  computational  Mach  number  and  wake  flow  angle.  Mach  number,  based 
on  stagnation  conditions,  was  redefined  with  respect  to  tangential  and  axial  velocity 
components. 


Program  plane.f 

^Hc5ie5|c*******5je5|c%*:je*jJe:Jc***5fc*3ie**5jc*5|c:icJic*5i«aJe5Je**si«*3k**5ie*****3ie%*******5le**5|« 

c  Modified  by  Lt.  William  Donovan  on  08  May  95  for  Master's  Thesis 
c  plane.f  reads  rvc3d  files  &  writes  ascii  files  for  plotting 
c  exit  plane  mach  number  and  flow  angles 
c  (2D  comparison  with  experimental  results) 
c  unit  1  =  input  xyz  file 
c  unit  3  =  input  q  file 
c  unit  4  =  output  residual  history 

c  xmit  7  =  output  (k  #1)  coimter,  theta,  pt/pO,Mach,  Angle 

c  unit  8  =  ouqjut  (k  #2)  counter,  theta,  pt/pO,Mach,  Angle 

c  xmit  9  =  output  (k  #3)  coxmter,  theta,  pt/pO,Mach,  Angle 

c  unit  10  =  output  (k  #4)  counter,  theta,  pt(pO,Mach,  Angle 

c  imit  11  =  ou^ut  (k  #5)  coxmter,  theta,  pt/pO,Mach,  Angle 

c 

parameter(ni=150,nj=3 1  ,nk=65) 
integer  um 

real  x(ni,nj,nk),y(ni,nj,nk),z(ni,nj,nk),  gama 

real  qq(5jii,njjik),resd(5000,5),ql,q2,q3,q4,q5,m(5) 

dimension  kk(5),ang(5),theta(5),tt0(5),vtnd(5),vmd(5), 

1  vznd(5),ptp0(5),deg(5) 
c 

c  k-values  are  hard-wired  below  (80.9%,  89.8%,  94.9%,  97.6%,  99.0%) 
c  (use  span.f  to  obtain  %-k  equivalency) 

C  (k  values  of  61  and  64  obtained  from  a  second  program  run) 
c 

data  kk/41,45,49,53,57/ 
c  read  grid  coordinates 
read(l  ,*)im,jm,km 

read(l,*)(((x(i,j,k),i=l,im),j=l,jm),k=l,km), 

1  (((y(i,j,k),i=l  ,im),j=l  jm),k=l  ,km), 

2  (((z(i,j,k),i=l  ,im),j=l  ,jm),k=l  ,km) 

^3{c:{e3k%*%3ie%:fc:je?fc:fe»{c^?ie:ic^?ie2ic9ic!ic3ie:je:jc:ic)ie^^sjc:jc;jcpje:jc:{e:{c3je:ic4e:|e^:|c:je9{c9{c?ic:k:||c){e%3ic:^:)c:je:{c9ic^9ie3ic:ic:ieHe 

c  read  restart  file 

read(3,*)imax,jmax,kmax 

read(3,*)fsmach,alpha,re,time 

c 

icheck=iabs(im-imax)+iabs(jm-jmax)+iabs(km-kmax) 

if(icheck.ne.O)then 

write(6,610)im,jm,km,imax,jmaxdcmax 

stop 

endif 

c 
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ooo  ooo  oooooo  ooo  ^ooo  oo 


read(3  *)((((qqG,id>k),i=l,im),j=l,jm),k=l,km),l=l,5) 


additional  residual  data 

read(3,*)itl,iil,phdeg,gama,om,nres,dum,dum,dum,dum 
read(3,*)((resd(nr,l),nr=  1  ,nres),l=l  ,5) 

Pressure  side  of  exit  wake  cut... 

=kk(3) 

i=l 

do  10  j=jmax,l,-l 
um=uni+l 
do  15  1=1,5 
k=kk(l) 

Rename  desired  q  values  for  use 

ql=qq(l,i,j,k) 

q2=qq(2,i,j,k) 

q3=qq(3,i,j,k) 

q4=qq(4,i,j,k) 

q5=qq(5,i,j,k) 

Compute  pl/pO  ratio 

ptp0(l)=gama*(gama-l)*q5+(gama-gama**2/2)/ql*(q2**2+q3**2+q4**2) 

Compute  T/Tt  ratio  and  mach  number 

tt0(l)=gama*(gama-l)/ql*(q5-0.5*(q2**2+q3**2+q4**2)/ql) 

m(l)=((l/tt0(l)-l)*2.0/(gama-l))**0.5 

Compute  Theta  and  Non-D  velocities  theta,  r  and  z 

theta(l)=asin(y(i,j,k)/(z(i,j,k)**2+y(i,j,k)**2)**0.5) 

vtnd(l)=-q4*sin(theta(l))+q3*cos(theta(l)) 

vmd(l)=q4*cos(theta(l))+q3*sin(theta(l)) 

vznd(l)=q2 

m(l)=(sqrt(vtnd(l)**2+vznd(l)**2))/ql 

dega)=theta(l)*57,296 

Compute  flow  angle  (degrees)  using  vmd  and  vznd 

ang(l)=atan2(vtnd(l),vznd(l))*57.296 
15  continue 

write(7,300)um,deg(l),ptp0(  1  ),m(  1 )  ,ang(l) 
write(8,300)um,deg(2),pt^(2),m(2),ang(2) 
write(9,300)um,deg(3),p^(3),m(3),ang(3) 
write(10,300)um,deg(4),ptp0(4),m(4),ang(4) 
write(l  1 ,300)um,deg(5),ptp0(5),m(5),ang(5) 

10  continue 


Repeat  above  on  other  side  of  wake  cut  at  imax 
c  j=l  not  used  (duplicate  location  across  wake) 
c  jmax  not  repeated  (same  location  as  j=31  at  i=l) 
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jmaxm=jmax-l 

i=im 


do  30  j=2,jmaxm,l 
ura=um+l 
do  25  1=1,5 
k=kk(l) 
ql=qq(l,i,j,k) 
q2=qq(2,i,j,k) 
q3=qq(3,i,jjc) 
q4=qq(4,i,j,k) 
q5=qq(5,i,j,k) 

ptp0(l)=gama*(gama-l)*q5+(gama-gama**2/2)/ql*(q2**2+q3**2+q4**2) 

tt0(l)=gama*(gama-l)/ql*(q5-0.5*(q2**2+q3**2+q4**2)/ql) 

m(l)=((l/tt0(l)-l)*2/(gama-l))**0.5 

theta(l)=asin(y(i,j,k)/(z(i,j,k)**2+y(i,j,k)**2)**0.5) 

vtndG)=-q4*sin(theta(l))+q3*cos(theta(l)) 

vmd(l)=q4*cos(thetaG))+q3*sin(theta(l)) 

vznd(l)=q2 

m(l)=(sqrt(vtnd(l)*  *2+vznd(l)  **2))/q  1 
dega)=theta(l)*57.296 
ang(l)=atan2(vtnd(l),vznd(l))  *57.296 
25  continue 

write(7,300)um,deg(  1 )  ,ptp0(  1  ),m(  l),ang(l ) 
write(8,300)um,deg(2),p^(2),m(2),ang(2) 
write(9,300)um,deg(3),p^0(3),m(3),ang(3) 
write(10,300)um,deg(4),ptp0(4),m(4),ang(4) 
writeO  1 ,300)um,deg(5),p^0(5),m(5),ang(5) 

30  continue 


c  residual  history  output  to  unit  4 

3|c  ^  Sfc  sjc  9|c  5jc  5jc  5jc  9|c  ^  9^  9{c  )|c  3|c  9|c  9|c  ^  9}c  jjc  9|c  9jc  9fc  Pjt  9|c  9|c  pjc  9|c  9^  Sjc  ^  djc  9{c  3|c  ^  9^  9}c  9jc  ^  9|c  9|C  ^  9fc  9|c  9{C  ^  9|C  ^  9jc  ^  ^  9jc  9^  9^  9)c  9^  9{c  9^ 

write(4,310)  l,(resd(l,l),l=l,5) 
do  40  j=2,nres 
it=10*0-l) 

*4p*^**|?;^*??/J^^^S’]y**?*********************************^ 

300  format(i5,5f8.3) 

310  format(i5,5(lx,el0.3)) 

610  formate  *****  warning  *****',/, 

1  ’  im,  jm,  km,  read  from  input', 3i5,’  do  not  match’,/, 

2  '  im,  jm,  km,  read  from  restart  file',3i5) 
stop 

end 
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Mach  Number  S  Mach  Number 


APPENDIX  H.  MACH  NUMBER  AND  FLOW  ANGLE  COMPARISONS 


Mach  Number  Comparison 

1.78  mm  Depth  (89.8%  Span) 


Computational  B  Experimental 


j  HI.  Mach  Number  Comparison  At  1.78  mm  Depth 


Mach  Number  Comparison 
0.42  mm  Depth  (97.6%  Span) 


Figure  H2.  Mach  Number  Comparison  At  0.42  mm  Depth 
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Flow  Angle  (Degrees)  |  ^  |  Mach  Number 


Wake  Flow  Angle  Comparison 

1.78  mm  Depth  (89.8%  Span) 


9-8-7-6-5-4-3-2  -1  0  1  2  3  4  5  6  7  8  9 


Wake  Position  (Degrees) 


Figure  H4.  Flow  Angle  Comparison  At  1.78  mm  Depth 
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Log  of  Residuals  S  Log  of  Residuals 


APPENDIX  L  CONVERGENCE  HISTORY 


RVC3D  Residuals  [0.5070  Pressure  Ratio] 


11.  0.5070  Pressure  Ratio  Convergence  History 


RVC3D  Residuals  [0.6041  Pressure  Ratio] 


Figure  12.  0.6041  Pressure  Ratio  Convergence  History 
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Log  of  Residuals  ^  Log  of  Residuals 


RVC3D  Residuals  [0.6815  Pressure  Ratio] 


13.  0.6815  Pressure  Ratio  Convergence  History 


RVC3D  Residuals  [0.8077  Pressure  Ratio] 


Figure  14.  0.8077  Pressure  Ratio  Convergence  History 
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Log  of  Residuals 


RVC3D  Residuals  [0.9054  Pressure  Ratio] 


Figure  15.  0.9054  Pressure  Ratio  Convergence  History 
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APPENDIX  J.  MULTIPLE  GRID  AND  SOLUTION  PLOTTING 


Program  mgrid.f  read  an  x,  y,  z  formatted  grid  file  (fort.l)  and  converted  it  into  a 
multiple  grid  file  format  (fort.21)  for  use  with  PLOTS D. 


Program  mgrid.f 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c  Program  to  read  a  single  PLOT3D  file  and  convert 
c  it  to  a  mgrid  file  format 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
real  x(250,65,65),y(250,65,65),z(250,65,65) 
real  yy(250,65,65),zz(250,65,65) 
read(l,*)idim,jdim,kdim 

read(l,*)(((x(i,j,k),i=l,idim),j=l,jdim),k=l,kdim), 

#  (((y(i,j  Jc),i=  1  ,idim),j=  1  ,jdim),k=  1  ,kdim), 

#  (((z(i  j  ,k),i=  1  ,idim),j=l  jdim),k=  1  ,kdim) 
write(*,*)’Done  reading  fort.!' 

c 

c  calculate  pitch 
c 

jm=jdim 

im=idim 

r2=y(  1  ,jm,  l)**2+z(  1  ,jm,  1  )**2 
sinp=(y(l,jm,l)*z(imom,l)-y(im,jm,l)*z(l,jm,l))/r2 
pitch=asin(sinp) 
c 

do  i=l,idim 
do  j=l,jdim 
do  k=ldcdim 

r=(z(i,j,k)**2+y(i,j,k)**2)**0.5 

theta=asin(y(i,j,k)/r) 

theta=theta+pitch 

y  y  (i  j  ,k)=r*  sin(theta) 

zz(i  ,j  ,k)=r*cos(theta) 

enddo 

enddo 

enddo 

c 

write(21)2 

write(2  l)(idim  jdim,kdim,igrid=l  ,2) 
c  do  10  igrid=l,2 

write(2  l)(((x(i,j,k),i=l  ,idim),j=  1  ,jdim),k=l  ,kdim), 

#  (((y(i  j  .k),i=  1  ,idim),j=  1  jdim)  ,k=  1  ,kdim), 

#  (((z(i,j  ,k),i=  1  ,idim),j= 1  jdim)  ,k=  1  ,kdim) 
c  10  continue 

write(2 1  )(((x(i  j  ,k)  ,i=  1  jdim)  j=  1  jdim),k=  1  ,kdim), 

#  (((yy(i.j.k)  j= 1  jdim)  j=  1  jdim)  jc=  1  .kdirn), 

#  (((zz(i  j  ,k)  j= 1  jdim)  j=  1  jdim)  jc=  1  ,kdim) 
stop 


Program  mq.f  read  a  solution  file  (fort.3)  and  converted  it  into  a  multiple  solution  file 
format  (fort.23)  for  use  with  PLOT3D. 


Program  mq.f 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c  Program  to  read  a  single  PLOT3D  file  and  convert 
c  it  to  a  Q_mgrid  file  format 
c 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
real  qq(250,65,65,5) 
read(3,*)idim,jdim,kdim 
read(3,*)fsmach,alpha,re,time 

read(3,*)((((qq(i,j,k,nx),i=  1  ,idim),j=  1  ,jdim),k=  1  ,kdim),nx=  1 ,5) 
c 

write(23)2 

write(23)(idim,jdim,kdim,igrid=l,2) 
do  10  igrid=l,2 
write(23)fsmach, alpha, re, time 

write(23)((((qq(i,j,k,nx),i= 1  ,idim),j=  1  ,jdim),k=  1  ,kdim),nx=  1 ,5) 

10  continue 
stop 
end 
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